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ABSTRACT 

We present a new family of spherically symmetric models for the luminous compo- 
nents of elliptical and spiral galaxies and their dark matter haloes. Our starting point 
is a general expression for the logarithmic slope a(r) — dlogp/dlogr from which most 
of the cuspy models yet available in literature may be derived. We then dedicate our 
attention to a particular set of models whose logarithmic slope is a power law function 
of the radius r investigating in detail their dynamics assuming isotropy in the velocity 
space. While the basic properties (such as the density profile and the gravitational 
potential) may be expressed analytically both the distribution function and the ob- 
servable quantities (surface brightness and line of sight velocity dispersion) have to be 
evaluated numerically. We also consider the extension to anisotropic models trying two 
different parameterization. Since the model recently proposed by Navarro et al. (2004) 
as the best fit to their sample of numerically simulated haloes belongs to the family 
we present here, analytical approximations are given for the most useful quantities. 

Key words: galaxies: kinematics and dynamics - galaxies: structure - stellar dy- 
namics - dark matter 



1 INTRODUCTION 

Investigating the kinematics and the dynamics of galactic 
systems is one of the classical task of astronomy. To com- 
plicate the issue, a galaxy is often made out of different 
subsystems and each one of them has its peculiarities that 
renders modeling the galaxy a quite difficult task. Along the 
years, this fascinating challenge has stimulated a lot of as- 
trophysicists to propose a plethora of different models that 
all share the common feature to be able to reproduce a set 
of observables suited for the particular kind of system one 
is interested to. 

The surface brightness is the easiest measurable prop- 
erty of a galaxy and is thus the first evidence a model builder 
has to take into account. This explains the large attention 
devoted to the search for models that, when projected on 
the plane of the sky, were able to reproduce the observed 
surface brightness, possibly mimicking the successful r 1//4 
law (de Vaucouleurs 1948). Illuminating examples of these 
kind are the models proposed by Jaffe (1983) and Hern- 
quist (1990). It is worth noting that analytical simplicity 
is another desirable property of a model. Hence, while it is 
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possible to deproject both the r 1 ^ 4 law and its generaliza- 
tion, the Sersic r 1 ^" profile (Sersic 1968), this gives rise to 
a model (Mazure & Capelato 2002) that is too complicated 
to be indeed useful. 

On the other hand, the dynamics of a galaxy may hardly 
be explained by its luminous components only thus claim- 
ing for the introduction of a dark halo whose shape, nature 
and structure are difficult to determine in a unique way. 
The impressive development in processing speed of modern 
computers have opened the way to more and more detailed 
simulations of structures formation that turned out to be a 
valid help in finding realistic and physically motivated mod- 
els for the dark halo. Actually, also this way is not free of 
problems mainly because of the unknown effect of numerical 
convergence and resolution on the asymptotic behaviour of 
the density profile. As a result, different simulations do not 
always give the same results and, moreover, it is not clear 
how these latter depend on the software used both to sim- 
ulate the system and to identify bound structures. While 
there is a general consensus that relaxed systems exhibit a 
density profile that is well described by a double power law 
with outer asymptotic slope —3, there is still an open con- 
troversy about the value of the inner asymptotic slope with 
proposed values mainly in the range ~ 1.0 — 1.5 (Navarro 
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et al. 1997; Tormen et al. 1997; Moore et al. 1998; Jing & 
Suto 2002; Power et al. 2003; Navarro et al. 2004). On the 
other hand, a similar controversy has arisen over the ques- 
tion whether such cusps are indeed observed in galaxies (see, 
e.g., Simon et al. 2003 and references therein). However, on 
galaxy scale, the effect of baryonic collapse and astrophys- 
ical feedback processes (such as supernova explosions) may 
alter significantly the dark halo structure thus biasing in a 
complicated way the interpretation of the observations. 

The great variety of models that have been proposed 
could disorientate giving the impression of a disordered col- 
lection of density profiles that are so different from each 
other that is somewhat surprising that all of them are able 
to describe the same kind of system. It is thus highly de- 
sirable to look for common features of all these models in 
order to see whether they are indeed different and indepen- 
dent descriptions of the same kind of system or rather they 
may be seen as particular cases of a more general family. 
A remarkable step in this direction is represented by the 
(a, f3, 7) - models proposed by Zhao (1996, 1997). Assuming 

spherical symmetry^, the mass density of the Zhao models 



p{r) = po 1 



l/a 



-1 -a(/3- 7 ) 



with r s a scale radius, ft and 7 the outer and the inner 
slope of the density profile respectively and a determining 
the width of the transition region. It is indeed possible to 
show that most of the previously proposed galaxy models 
(both cuspy and cored) are obtained by suitably setting the 
three parameters (a,/3,7). Moreover, the Zhao models also 
forms a complete set for constructing general galaxy models 
or solving Poisson equation in the non spherical case. 

There is still another possible approach to smooth out 
the differences among spherically symmetric cuspy models. 
Qualitatively, we can characterize a model by its asymptotic 
behaviours for r — » and for r — > 00. A straightforward way 
to do this is to assign the logarithmic slope of the density 
profile, i.e. a(r) = dlog p/d log r. While the usual approach 
consists in giving the mass density and then evaluating a(r), 
the inverse way is also possible. We first choose a general 
expression for the logarithmic slope and then evaluate the 
corresponding density profile by a simple integration. As 
we will show later, this method allows to reduce most cupy 
models to some few classes on the basis of the shape of 
a(r). Moreover, it is worth stressing that a(r) may be easily 
estimated by numerical simulations of structure formation 
and is less affected by systematic effects difficult to control. 

The paper is organized as follows. In Sect. 2 we intro- 
duce a general parameterization of the logarithmic slope 
a(r) and delineates four classes of models that are worth 
to be investigated. In particular, we present in this paper 
a detailed study of the kinematics and dynamics of models 



T Without loss of generality, we will only consider spherically 
symmetric models. Although real galaxies are seldom spherically 
symmetric, they are indeed axisymmetric, while dark haloes are 
usually assumed to be spherically symmetric being their true 
shape unknown. Moreover, introducing an ellipticity in the den- 
sity profile complicates evaluating the model properties without 
altering the main results. 



with power law logarithmic slope. Sect. 3 is devoted to the 
estimate of the basic properties : we evaluate the mass den- 
sity, mass profile, circular velocity, gravitational potential 
and isotropic velocity dispersion of this class of models. The 
dynamical properties are fully characterized by the distri- 
bution function and the density of states that are estimated 
in Sect. 4 assuming isotropy in the velocity space. In Sect. 5, 
we evaluate the observable quantities, namely the surface 
density and luminosity weighted projected velocity disper- 
sion, while, in Sect. 6, we give off the hypothesis of isotropy 
and investigate the effect of anisotropy on the radial velocity 
dispersion and the distribution function. We summarize and 
conclude in Sect. 7. 



2 THE LOGARITHMIC SLOPE 

The usual approach to modeling galactic systems (both the 
luminous components and the dark haloes) starts from the 
mass density profile and proceeds by evaluating the main dy- 
namical quantities. The increasing resort to numerical sim- 
ulations of structure formations has however drawn more 
attention to the logarithmic slope of the radial profile de- 
fined as : 



a(r) = 



dlogp 
dlog r 



(1) 



It is indeed easier to compare different galactic models 
based on their corresponding logarithmic slope. Moreover, 
the knowledge of a(r) renders quite immediate to study the 
asymptotic behaviours (towards the centre and/or the infin- 
ity) of the density profile. Furthermore, it seems that this 
quantity is better constrained than the density profile by 
numerical simulations. 

Motivated by these considerations, it is interesting to 
look for an expression of a(r) that is as more general as 
possible in order to build spherical galactic models that may 
be grouped under the same class by mean of the property 
that they share the same logarithmic slope. To this aim, a 
useful proposal for a(r) is : 



l+a(r/r a )+b(r/r s f 
' 1 + c(r/r s ) + d(r/r.) a 



(2) 



with r s an arbitrary scaling radius. It is possible to show 
that most of the galaxy models proposed up to now in lit- 
erature may be obtained by inserting Eq.(2) into Eq.(l) 
and integrating with respect to r with the boundary con- 
dition p(r s ) = p s after having suitably set the parameters^ 
(/?, a, b, c, d, 7). Unfortunately, the result we get by letting 
completely unspecified the above parameters turns out to be 
extremely complicated involving the exponential of a com- 
bination of hypergeometric and trigonometric functions. As 
a consequence, it is not surprising that also the simplest 
dynamical quantities (such as the mass profile) are analyti- 
cally impossible to derive so that Eq.(2) is useless for astro- 
physical applications. We have thus to look for general, but 
simpler expressions for the logarithmic slope a(r). 



t In Appendix A we give the values of the parameters 
(/3, a, b, c, d, 7) for some popular models. 
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To this aim, let us observe that most of the cuspy mod- 
els frequently used in literature may be deduced by integrat- 
ing Eq.(2) with b = d — so that, heron, we will restrict our 
attention to models having the following expression for the 
logarithmic slope : 



, , a 1 + a(r/r s ) ( r 

air) = -p x )-!-. — f x — 

V ' l + c(r r s ) \r s 



(3) 



Let us first consider the case c =fc 0. Without loss of gener- 
ality, we may set c = 1 thus obtaining : 



i \ a 1 + a ( r / r s) 
W 1 + (r/r 3 ) 



(4) 



If c / 1, we may rescale all the results obtained starting 
from Eq.(4) by replacing: 

r s — » r s /c , (3 — > P/c 1 , a — » a/c . 

For 7 = 0, Eq.(4) reduces to : 

1 + a(r/r s ) 



a{r) = -/)x 



1 + (r/r s ) 



(5) 



that may be considered as a generalization of a lot of double 
power law models (see Table 1 in Appendix A) . On the other 
hand, by setting a = c = and 7 / 0, we get : 



«(r) = -/*x(£)' 



(6) 



For (/3, 7) = (2,0.17), this reduces to the model recently 
proposed by Navarro et al. (2004, hereafter N04) on the 
basis of a set of high resolution numerical simulations of 
dark matter haloes. Finally, for c = and 7 / 0, we may 
set a — 1 and obtain a fourth class of models with : 



a{r) = -/)x (1 + r 



(7) 



that may be generalized to the case a 7^ 1 by the following 
replacements : 

r s — » r s /a , /? — > /3/a 7 . 

Summarizing, we have defined four set of models charac- 
terized by logarithmic slopes given by Eqs.(4), (5), (6) and 
(7) respectively. We have checked that, while it is possible 
to get an analytical expression (at least, in terms of spe- 
cial functions) for the density profile in all four cases, the 
mass profile is analytical only for models with a(r) given by 
Eq.(5) and Eq.(6). Our aim here is to study galaxy mod- 
els that are both general and analytically amenable so that 
they may be easily compared to observational data (such as 
those on the rotation curve). That is why only Eqs.(5) and 
Eq.(9) are worth to be considered in detail. Moreover, it is 
easy to show that integrating Eq.(5) gives rise to a class of 
models that is only a subset of the more general family of the 
Zhao models. This is partially true also for models with a(r) 
given by Eq.(6) that may indeed be obtained as a limiting 
case. However, the dynamical properties of these particular 
models have not been investigated in detail by Zhao. More- 
over, the result of Navarro et al. (2004) quoted above is a 
strong motivation to dedicate much attention to these mod- 
els which we will refer to in the following as PoLLS {Power 
Law Logarithmic Slope) models. 



3 BASIC PROPERTIES 

The logarithmic slope of the radial density profile for PoLLS 
models is given by Eq.(9) and is characterized by four pa- 
rameters, namely the scaling radius r s , the characteristic 
density p s and the two slope parameters (/?, 7). Actually, it 
is possible to reduce the number of parameters by redefin- 
ing the scale radius. To this aim, let us evaluate r-2, the 
radius at which the logarithmic slope equals the value of the 
isothermal sphere. By solving a(r_ 2 ) = —2, we easily get : 



r-2 



1/7 



Replacing r s with r_2, Eq.(9) may be written as: 



a(r) = -2 

r-2 



(8) 



(9) 



so that the parameter f) may indeed be eliminated. It is 
interesting to look at the asymptotic behaviours of a(r) : 



lim a(r) 

r— »0 



, 



lim a(r) = —00 . 

r— >oo 

The logarithmic slope does not diverge in the centre so that 
the model is not singular which is a nice feature, while a 
is monotonically decreasing thus suggesting an exponential - 
like decrease of the mass density. Actually, the density profile 
turns out to be : 



p(r) 



P-2 exp 



(10) 



so that p remains finite for r — > 0, which is what usually 
happens in the case of cored models as, e.g., the non singu- 
lar isothermal sphere (Binney & Tremaine 1987; Hinshaw & 
Krauss 1987). The exponential decrease of p is coherent with 
the asymptotic limit of a for r — > 00. As a general remark, 
note that PoLLS models, loosely speaking, may be consid- 
ered as a one parameter family with 7 as ordering parameter 
and (r_2,p— 2) scaling quantities. 

A strong constraint on 7 comes from the evaluation of 
the mass profile : 



M(r) = 47r / r' 2 p(r')dr' = 4np-2r s 



f 

Jo 



exp 



7 



7-1-, 



having defined x = r /r-2- This integral is defined only for 
7 > giving : 

r(3/ 7 ) - r(3/ 7 , 2^/7) 



M(x) = Mtot x 



r(3/ 7 ) 

where Mtot is the total mass given by : 

3/7 

r(3/ 7 ) exp (2/7) 



4np. 



- 2 r 3 _ 2 



7 



(11) 



(12) 



and r(a;) and F(a, x) are the usual F function and the in- 
complete r function respectively. In Fig. 1, we report the 
logarithm of M n = M(x)/M to t for three different choices of 
7. Except in the very inner regions (log a; < — 1), for a fixed 
value of x, the higher is 7, the higher is the ratio among the 
mass within x and the total one thus indicating that models 
with higher values of 7 converges more quickly to the total 
mass. It is worth noting, however, that M n converges to 1 
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Figure 1. Logarithm of the normalized mass as function of the 
logarithm of the scaled radius x for PoLLS models with three 
values of 7, namely 7 = 0.085 (short dashed), 7 = 0.17 (N04 
model, solid) and 7 = 0.34 (long dashed). 
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Figure 2. Logarithm of the scaled total mass M s t = 
Mtot/4 : np-2r 3 _2 as function of the logarithm of the slope 7. 



only slowly, i.e. for values of r » 5r_ 2 for the models we 
have considered. As a result, the scaled half mass radius Xhm, 
implicitly defined as M(ij, m ) / 'Mtot = 1/2, turns out to be 
quite large. For instance, for the N04 model it is x hm ~ 10, 
while for 7 6 (0.06, 1.06) Xhm is well approximated (within 
5%) as : 



\ogx hm ~ 0.1467" 



(13) 



as we have checked numerically. As regard the total mass, 
this is decreasing function of the slope parameter 7 as it 
is shown in Fig. 2. Qualitatively, this can be explained by 
noting that models with higher 7 are more concentrated (i.e. 
the density decreases faster) so that less mass is contained 
within a given volume with the result that Mtot is smaller. 

Being the model spherically symmetric, we do not need 
to explicitly compute the gravitational potential to evaluate 
the rotation curve. It is simply : 



Voir) = y/GM(r)/r 
that, in our case, reduces to : 
V c (x) — V-2 x V(x) 
having defined : 



(14) 
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Figure 3. The normalized circular velocity v n (x) = 

v ci%)/ \J GMtot/r-2 as function of the logarithm of the scaled ra- 
dius x for PoLLS models with three values of 7, namely 7 = 0.085 
(short dashed), 7 = 0.17 (N04 model, solid) and 7 = 0.34 (long 
dashed). 



v-2 = v c (r = r-2) 



r-2 



r(3/ 7 ) - r(3/ 7 , 2/7) 



r(3/ 7 ) 



Vix) 



1 r(3/ 7 ) - r(3/ 7 , 2^/7) 
x r(3/ 7 )-r(3/7,2/ 7 ) 



,(15) 



(16) 



Looking at Fig. 3, where the normalized normalized velocity 
v„(x) = v c (x)/W GMtot I r-2 is shown, we see that higher 
values of 7 gives rise to higher peak velocities. The position 
x m ax of the peak is a function of 7 and is well approximated 
(within 5%) as : 



I.68O7" 



(17) 



For the N04 model, it is x max — 2.1, while models with 
lower 7 peaks farther from the galaxy centre. The rota- 
tion curve is asymptotically keplerian. Indeed, for x — > 00, 
r(3/7, 2x J /~/) — > so that the fractionary term in Eq.(16) 
goes to unity and the keplerian term 1/x dominates. Actu- 
ally, this is an obvious consequence of the finite total mass 
of the model. Note, however, that the keplerian behaviour is 
attained only for large values of x so that v c is quite slowly 
declining and, within the observational errors, may be also 
compatible with a flat rotation curve. 

Let us now derive the gravitational potential for PoLLS 
models. To this aim, the Poisson equation : 

V 2 $ = 47rGp(r) 

has to be solved. For spherically symmetric system, the gen- 
eral solution is (Binney & Tremaine 1987) : 



$( r ) = „^£W _ 4ttG^ pir')r'dr' 



(18) 



Using Eq.(ll) for the mass profile and Eq.(10) for the den- 
sity law, we get : 



$(x) = -$_ 2 (p-2,r_2,7) x 



with : 

GMtot % 

$-2 = x T 1;7 , 

r-2 



H x n) 
^(i;7) 



(19) 



(20) 
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Figure 4. The scaled central potential <3?o s as function of the 
logarithm of the slope parameter 7. We have defined $os = 
-*(0)/(GM fo( /r_ 2 ). 
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Figure 5. The velocity dispersion cr (in km/s) as function of 
the radius r (in kpc) for PoLLS models with three values of 7, 
namely 7 = 0.085 (short dashed), 7 = 0.17 (N04 model, solid) and 
7 = 0.34 (long dashed). The fiducial values M to t = 5 X 10 11 Mq 
and r_2 = 10 kpc have been used. 



+ 



r(3/ 7 )-r(3/ 7 , 2x^/7) 



xr(3/ 7 ) 

r(2/ 7 ) + r(2/ 7 ,2^/ 7 ) 
r(3/ 7 ) 



1/7 



(21) 



Being the total mass finite, Eq.(18) easily shows that 
the gravitational potential vanishes at infinity. However, a 
caveat is in order here. If we consider the limit for x — > 00 of 
Eq.(19), we find $(00) = GM tot /r- 2 x r(2/ 7 )/r(3/7). Since 
the potential is defined by Eq.(18) up to an arbitrary addi- 
tive constant, we could add this value to Eq.(19) in order 
to have <&(x) vanishing at infinity as it is indeed the case. 
At the opposite end, having the model a finite force at the 
centre, the gravitational potential does not diverge in the 
origin. By applying the limit for x — > in Eq.(19), we get : 



*(0) 



GMto 
r-2 



1/7 



r(2/ 7 ) 
r(3/ 7 ) 



(22) 



which is plotted in Fig. 4. The central potential turns out 
to be an increasing^ function of the slope parameter 7. As 
a consequence, the depth of the potential well is greater for 
larger values of 7 and, for these models, particles are more 
attracted toward the centre. Therefore, for a given x, the 
gravitational force and hence the number of stars within a 
distance x from the centre is larger for larger values of 7 
thus explaining why M (x) / M to t is an increasing function of 
the slope parameter. 

Another dynamically interesting quantity that may be 
straightforwardly evaluated is the velocity dispersion. As- 
suming isotropy in the velocity space, this is given by (Bin- 
ney & Tremaine 1987) : 



1 



P(r) J r 

Inserting Eqs.(10) and (11) into the above general expres- 
sion, we get : 



GM tot exp (-2/7) 
T(3/ 7 ) 



X I a {x;i) 



with 



J X 



rr(3/ 7 )-r(3/7,2T/7)]e 



-(2/ 7 ) S ^-2^ 



(23) 



(24) 



This quantity may not be evaluated analytically, but it is 

straightforward to estimate numerically^ for fixed values 
of the slope parameter 7. Fig. 5 shows the velocity disper- 
sion for three values of 7 using M to t = 5 x 10 11 M and 
r_2 = 10 kpc. The result may be simply scaled to other 
values of (M tot ,r_2) by noting that, because of Eq.(23), 
a oc Mtot /r-2- As yet observed for the rotational velocity, 
for a given x, a is higher for higher values of 7. In both cases, 
we have thus a degeneracy between the total mass and the 
slope parameter since both these quantities increase v c and 
cr. However, the degeneracy may be broken by determining 
where the rotation curve peaks since this only depends on 7. 
As a further remark, we note that a decreases in the outer 
regions of the system thus mimicking well what has been re- 
cently observed in intermediate luminosity elliptical galaxies 
(Romanovsky et al. 2003) . 



4 THE DISTRIBUTION FUNCTION 

All the kinematical information on a collisionless system 
may be derived by the phase space distribution function 
(hereafter DF) /(r,v) that, for a non rotating spherically 
symmetric isotropic system, depends only on the binding 
energy for unit mass. Following Binney & Tremaine (1987), 
we define : 



*(r) = -$(r) + $0 , 8{r) 



§ Note that Fig. 4 report — <I>(0) so that a monotonically decreas- 
ing curve means an increasing central potential. 



» In Appendix B we give analytical expressions that approximate 
very well some of the dynamical quantities we compute in this 
paper for the special case of the N04 model. 
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Figure 6. The logarithm of the reduced DF / rec j as function of 
the logarithm of the dimcnsionlcss binding energy £ for PoLLS 
models with three values of 7, namely 7 = 0.085 (short dashed), 
7 = 0.17 (N04 model, solid) and 7 = 0.34 (long dashed). 



Figure 7. The logarithm of the reduced density of states g re d as 
function of the logarithm of the dimensionless binding energy £ 
for PoLLS models with three values of 7, namely 7 = 0.085 (short 
dashed), 7 = 0.17 (N04 model, solid) and 7 = 0.34 (long dashed). 



and set $0 = so that the DF vanishes for £ < 0. The DF 
may then be obtained by means of the Eddington formula 
(Binney & Tremaine 1987) : 



/(£) = 



1 



d 



dp d* 



y87r 2 d£ J d* v/£^* 



In principle, one should express p as function of * by first in- 
verting the relation ^ = 't'(x). Unfortunately, Eq.(19) may 
not be inverted analytically so that we have to resort to a dif- 
ferent approach changing variable from $ to x in the above 
Eddington formula. Using Eq.(10) and (19), we rewrite the 
Eddington formula as : 



f{£) 



2p_ 



r- 1 exp[-(2/7)(r-l) 



<(25) 



7r2$_ 2 V8*-2 d£ J xi s) y/£ - 

having introduced the dimensionless quantities : 

£ = £/®-2 , * = */$_ 2 . 

In Eq.(25) the lower limit of integration is implicitly defined 
by the equation : 



*[x(£)] 



(26) 



that is easy to solve numerically. The dimensionless energy 
£ ranges between |$(oo)/$_2|/ and |$(0)/4>-2| that, in our 
case, reduces to < £ < |$(0)/$_2| with $(0) given by 
Eq.(22). For values of £ in this range, we may numerically 
determine the DF of the PoLLS models by first setting the 
value of the slope parameter 7, then solving Eq.(26) with re- 
spect to x(£) and finally using Eq.(25). A caveat is in order 
here. For £ running in its range, x changes between and 
00. Actually, when numerically solving Eq.(26), there are 
technical problems approaching the lower and upper ends of 
the range for the dimensionless binding energy. As a conse- 
quence, we compute the DF only over a limited range. The 
result is shown in Fig. 6 where we plot the logarithm of the 
reduced DF defined as : 



!r,i(£ ) = f{£) x 



2p-2 



2 $- 2 ^8¥~ 



It is worth noting that the reduced DF is well approximated 
by a single power law for a large range of values of the bind- 
ing energy £ . This is apparent from Fig. 6 where the loga- 



rithm of the reduced DF turns out to be a linear function of 
log £ for log £ . However, both the exponent of the power law 
(i.e., the slope of the linear part of the plot) and the range 
over which this approximation works sufficiently well depend 
in a complicated way on the parameter 7. In particular, the 
range over which the reduced DF is a power law reduces with 
increasing 7. As a general remark, we note that the shape 
of the logarithm of the reduced DF is the same for different 
values of 7 with a linear part followed by a steep increase. In 
particular, for a fixed value of £, the higher is 7, the larger 
is the reduced DF. It is important to stress, however, that 
the power law approximation does not hold at the low and 
high energy limits. Actually, for £ approaching |$(0)/$_2|, 
the DF diverges (as shown by the steep increase in Fig. 6) 
and the divergence is strongest for smaller 7. On the other 
hand, f re d drops to zero much steeper than a power law for 
£ — > (not shown in the plot), but, unfortunately, we have 
been unable to find an approximated formula relating the 
steepness of the decrease with the slope parameter 7. 

The DF is a useful quantity, but does not tell anything 
about the number of stars per unit binding energy. The 
quantity that gives this information is the differential en- 
ergy distribution M{£) that, for an isotropic and spherically 
symmetric system, may be written as Af(£) = f(£)g(£). 
Since we have yet determined the DF, computing the differ- 
ential energy distribution reduces to the evaluation of g(£ ), 
that is the density of states per unit binding energy given 
by (Binney & Tremaine 1987) : 



g(£) = 16\/27r 2 



p 00 


2 dr 


Is 


T d* 



VW~£ d* 



Since \& is a monotonic decreasing function of r, dr /d^ < 
everywhere and we can rewrite the above relation as : 



/•*(£) 

g{£) = 16 v /2$-27r 2 r 3 _2 / VW~£ f 
Jo 



(27) 



with x(£) the same as in Eq.(26). It is convenient to define 
a dimensionless reduced density of states as : 



g re d(£) 
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Figure 8. The logarithm of the reduced differential energy dis- 
tribution N re d as function of the logarithm of the dimension- 
less binding energy £ for PoLLS models with 7 = 0.085 (short 
dashed), 7 = 0.17 (N04 model, solid) and 7 = 0.34 (long dashed). 
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Figure 9. The logarithm of the scaled surface density E n = 
E/E_2 as function of the logarithm of the scaled radius y for 
PoLLS models with three values of 7, namely 7 = 0.085 (short 
dashed), 7 = 0.17 (N04 model, solid) and 7 = 0.34 (long dashed). 



As is shown in Fig. 7, the density of states may be approxi- 
mated as a power law over a large range of £ as yet observed 
for the reduced DF, but now the slope is negative coherently 
with what is expected. The dependence on 7 is similar to 
that of the reduced DF with the only difference that the 
reduced density of states is a decreasing rather than an in- 
creasing function of the binding energy. 

Finally, we show in Fig. 8 the differential energy distri- 
bution by plotting the logarithm of the reduced quantity : 



M red (£)=M(£) x 



l6p- 2 rL 2 



Af turns out to be a decreasing function of the binding en- 
ergy well approximated by a single power law (i.e. a linear 
fit in logarithmic space) over almost the full range of £. The 
abrupt cutoff at high binding energies may be qualitatively 
explained by noting that high values of £ are reached in the 
centre of the mass distribution where the density takes a 
finite value. As a result, the differential energy distribution 
may not go down forever, but has to be truncated. At the 
other extreme, since the mass density formally vanishes only 
at infinity, it is always possible to find stars at larger and 
larger distances with lower and lower binding energies so 
that M{£ ) may diverge for £ going to zero. This behaviour 
is apparently opposed to what is shown in Fig. 8 where a 
divergence appears only for the model with 7 = 0.085, while 
the Nred is approximately constant for 7 = 0.17 and de- 
creasing for 7 = 0.34. However, we have checked that these 
unexpected behaviours are the result of numerical problems 
originating for low values of the binding energy. 



5 OBSERVED PROPERTIES 

The density law in Eq.(10) is a generalization of the N04 
model that has been proposed to describe the dark haloes 
coming out from high resolution numerical simulations. Ac- 
tually, there is no a priori reason against using PoLLS mod- 
els to describe also the luminous components of galaxies (in 
particular, ellipticals and bulges of spiral galaxies). It is thus 
interesting to ask what is the surface mass density of these 
models since this may compared with the surface brightness. 



Moreover, this quantity directly enters the study of the lens- 
ing properties of a given model. It is thus worth computing 
the surface density of PoLLS models starting from the defi- 
nition : 



£(i?) 



t' R p(r)rdr 
~ Jo Vr 2 - R 2 



with R the distance from the centre in the equatorial plane 
of the system. By inserting Eq.(10) into the above relation, 
we get : 



E(y) = £-2 x 



5(1,7) 

with y = R/r-2 and we have defined : 
ry e-K V ,-{2/ 1 ){x 1 - 1 )xdx 



ry 

S(y,-y)= / 
Jo 



\A 2 - y 2 



S_ 2 = 2p_ 2 r_ 2 <S(l,7) . 



(28) 

(29) 
(30) 



<S(j/,7) may not be evaluated analytically, but it is straight- 
forward to compute numerically for fixed value of the slope 
parameter 7. The result is shown in Fig. 9 where we use a 
logarithmic scale to make easier the visualization. It is worth 
noting that the slope parameter 7 plays a key role in deter- 
mining the behaviour of E in the inner regions, while has 
a still significant but less dramatic effect for large values of 
y. In particular, we find that the higher is 7, the lower is 
the surface density for a given x. This is particularly evi- 
dent in the very inner regions as it is seen comparing the 
flat part of the curves in Fig. 9. It is worth stressing that 
the surface density may not mimic neither the usual r 1 / 4 
law nor its generalization r 1 ^™. This is not a serious draw- 
back of PoLLS models, however, since they were introduced 
to describe dark haloes rather than luminous components. 
It is nonetheless interesting to compute also the integrated 
surface density that is given as : 



f 

Jo 



S(y) = 2tt / T,(R')dR' = 27r£_ 2 r 2 



S(y',"f) 
5(1,7) 



dy (31) 



that is shown in Fig. 10. The dependence on 7 is further 
weakened by the integration process with the result that 
S(y) is almost independent on the slope parameter up to y ~ 
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Figure 10. The normalized integrated surface density S n = 
S/ (27rE_2»*^ 2 ) as function of the scaled radius y for PoLLS mod- 
els with three values of 7, namely 7 = 0.085 (short dashed), 
7 = 0.17 (N04 model, solid) and 7 = 0.34 (long dashed). 



2. Note, however, that both S n and S„ = S , /(27rE_2r 2 _2) are 
increasing function of 7 as it is expected since models with 
smaller values of 7 are more concentrated, i.e. the projected 
mass within a given distance from the centre is larger. 

In Sect. 3, we have evaluated the velocity dispersion, 
but indeed it is not this quantity that is measured. Actually, 
from the observations it is possible to estimate the velocity 
dispersion projected along the line of sight and weighted 
with the surface luminosity density. This is given by : 



2 

O"(os — 



I(R) 



J R 



pstars (r) 



GM T (r)Vr 2 - K 2 



dr 



(32) 



where I(R) is the total surface luminosity density, pstars 
the mass density of the stellar component and Mr the total 
(stellar plus dark matter) mass profile. To be more realistic, 
we average over a finite radial bin (i?i,7?2) as follows (see, 
e.g., the appendix in Dalai & Keeton 2003) : 



<?bv 



with : 



A{x) = 



F(x) 



J£ 2 Pstars(r)GM T (r) [F(R 2 /r) - F(Ri/r)] dr 
2 P stars(r) [A(R 2 /r) - A(Jii/r)] rdr 

arcsin(^) x < 1 

7r/2 X > 1 



xy/1 — x 2 + arcsin(:r) — ty/2 x < 1 



(33) 



(34) 



(35) 







x > 1 



A caveat is in order here. In principle, I(R) is not equal to 
T,(R). Actually, we have to distinguish among three possi- 
bilities. First, let us suppose that we use a PoLLS model 
for the dark halo of a galaxy in which case E(-R) does not 
enter in I{R), while the mass profile of the model still en- 
ters in determining the velocity dispersion through Mt in 
Eqs.(32) and (33). On the other hand, a PoLLS model may 
also be used to fit the surface brightness profile of an ellip- 
tical galaxy so that we need both E, p and M to evaluate 
ai s and a bin and also choose a model for the dark halo. 
Finally, if the galaxy is a spiral one, we have a three com- 
ponent system thus complicating a lot describing the whole 



Figure 11. The projected luminosity weighted velocity disper- 
sion averaged over a radial bin cr oi „ (in km/s) as function of 
the distance R (in kpc) for PoLLS models with three values of 7, 
namely 7 = 0.085 (short dashed), 7 = 0.17 (N04 model, solid) and 
7 = 0.34 (long dashed). The fiducial values M tot = 5 X 10 11 M Q 
and r_2 = 10 kpc have been used. 



system and computing the velocity dispersion. In order to 
not introduce degeneracy among the model parameters and 
to better investigate the properties of the PoLLS models, we 
restrict our attention to elliptical galaxies and assume that 
it is possible to use a PoLLS model to fit the surface bright- 
ness. Moreover, we arbitrarily set R2 = l.lRi which implies 
that the larger is the distance from the centre, the larger is 
the bin width thus allowing to reproduce qualitatively the 
need of averaging over larger bins in the outer regions to 
ameliorate the signal to noise ratio. 

Under these assumptions, pstars is given by Eq.(10) and 
M T (r) by Eq.(ll) so that Eq.(33) becomes: 

VbiniVl) = G ^° l . r X hin(yi) (36) 



2r_ 2 F(3/ 7 ) 



with : 



IbiniVl) = A4m(2/l)/£>6in(2/l) , 

f V1 [T(3/ 7 ) - r(3/ 7 , 2x77)] * 



(37) 



J VI 



exp[-(2/7)(^-l)] x 
[F(l.ly 1 /x)-F(y 1 /x)]dx , 

a;exp[-(2/ 7 )(a; 7 - 1)] x 

L4(1.1j/i/ie) - A(yi/x)\ dx . 



(38) 



(39) 



Not surprisingly, we are unable to give an analytical ex- 
pression for I bin , but it is not difficult to evaluate it nu- 
merically for fixed value of the slope parameter 7. To get 
some interesting examples, we use the same fiducial model 
adopted for computing a in Fig. 5 obtaining the plot in 
Fig. 11. Noting that, by virtue of Eq.(36), atin is propor- 
tional to Mtot/ r -2, it is immediate to scale the results for 
the fiducial model to other values of these two parameters. 
It is worth noting that there is not a clear trend of aun with 
7 in the very inner regions (j/i < 0.05 not visible in the plot) 
as a result of the complicated interplay between the density 
and the mass profiles due to the fact that M{r) enters only 
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in Mbin{yi)- However, for larger values of yi, we recover the 
same behaviour yet observed for the velocity dispersion a 
shown in Fig. 5 with aun increasing with 7. Fig. 11 does not 
extend further than y — 0.1. While there are no technical 
difficulties in computing aun for y > 0.1 thus obtaining a 
declining projected velocity dispersion, Eqs.(36) - (39) rigor- 
ously do not hold anymore in this regime since they have 
been obtained under the hypotheses that there is no dark 
matter (or that its contribution may be neglected) and a 
PoLLS model may describe the luminous component. Both 
these hypotheses (in particular, the former one) introduce 
systematics that quickly increase approaching the outer re- 
gions of the galaxy so that we have preferred to only evaluate 
abi n in the very interiors to reduce these errors. 



6 ANISOTROPIC MODELS 

The models in the family we have described have an isotropic 
dynamical structure, i.e. we have assumed that the velocity 
dispersion is the same along the three axes of the velocity 
ellipsoid. Actually, this hypothesis is hardly verified in real 
systems so that it is important to extend our models to take 
into account a possible anisotropy in the orbital structure. 
The radial velocity dispersion oy is determined by solving 
the Jeans equation that now writes : 



d(pal) „3 ff (r) , . j. . , GM(r) 

with 



_2 

1 _ £n 

1 



(40) 



(41) 



where at is the tangential velocity dispersion. In the 
isotropic case, at = a r = a and Eq.(23) holds. In the general 
case, the solution of Eq.(40) is (Binney & Tremaine 1987; 
Mamon & Lokas 2004b) : 



V(r')p(r') 



(42) 



(43) 



with r)(r) related to the anisotropy parameter /3 a as : 
dlogr 

The estimate of the anisotropy parameter j3 a is quite dif- 
ficult so that different choices are possible. A widely used 
parameterization is the Osipkov- Merrit one (Osipkov 1979; 
Merritt 1985) : 

„2 



3 a {x) = 



(44) 



ry2i I -y £ 

x 1 x OM 

that is characterized by the anisotropy radius Tom = 
t 3 xom- Solving Eq.(43) with /3 CT given by Eq.(44) and in- 
serting the result and Eqs.(10) and(ll) into the general ex 
prcssion (42), we get : 



2 GMtot 
a r — X 

r-2 



{r(3/ 7 )(x 2 + xIm) exp [-(2/ 7 )(^ - I)]}" 1 
x Iom{x;j,xom) (45) 



with : 

Iom - 



f 

■J X 



[T(3/7)-r(3/7,2r/7)] x 
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Figure 12. The radial velocity dispersion a (in km/s) of the N04 
model as function of the radius r (in kpc) for the isotropic model 
(solid), the Osipkov - Merritt models with xqm = 0-5 (short 
dashed) and xqm = 1 (long dashed) and the Mamon & Lokas 
model with xml = 0.18 (continuous dashed). The fiducial values 
Mtot = 5 X 10 11 Mq and r_2 = 10 kpc have been used. 



exp[-(2/ 7 )(r-l)] 



£ + x OM 



(46) 



Even if widely used, it has been recently shown that the 
Osipkov - Merritt parameterization does not fit well the 
anisotropy parameter measured in cosmological N-body 
simulations (Diaferio 1999; Colin et al. 2000; Rasia et 
al. 2004; Diemand et al. 2004). Indeed, Mamon & Lokas 
(2004a, b) have shown that the following parameterization 
works better : 



2{x + xml) 



(47) 



with xml = rML/r s and tml a characteristic anisotropy 
radius. With the Mamon & Lokas parameterization of the 
anisotropy parameter, the radial velocity dispersion turns 
out to be : 

GM to t 



a r {x) = 



r-2 



x {r(3/7)(x + x M l) exp [-(2/ 7 )(x 7 - l)]}" 1 



x Iml{x\~i,x M l) 



(48) 



with : 



-/ 

J X 



[r(3/ 7 )-r(3/7,2C/7)] x 



exp[-(2/ 7 )(r 



1)] i+|£££ di 



(49) 



Actually, it is worth noting that some of the simulations 
considered by Mamon & Lokas are tailored on cluster scale 
haloes (Rasia et al. 2004) or on the subhaloes (Diemand 
et al. 2004) so that it is at least premature to reject the 
Osipkov - Merritt model for galaxy scale systems. Hence, in 
Fig. 12, we report the radial velocity dispersion a r for the 
N04 model for the isotropic model, the Osipkov - Merrit and 
the Mamon & Lokas anisotropic parameterizations. As ex- 
pected, anisotropic models predict an higher velocity disper- 
sion with the Osipkov - Merritt scheme being more effective 
in raising a r . In particular, it turns out that the lower is 
xom, the higher is the the radial velocity dispersion that 
may be qualitatively explained by noting that the model be- 
comes more and more isotropic as xqm increases. We have 
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Figure 13. The logarithm of the reduced DF of the N04 model as 
function of log (with © = Q) for three values of the dimension- 
less anisotropy radius xqm > namely xqm = 0-5 (short dashed), 
%OM = 1-0 (solid) and xqm = 2.0 (long dashed). 



checked that these results do not depend on 7. Note that 
there is a degeneracy between mass and anisotropy so that 
it is difficult to choose among isotropic, Osipkov - Merritt 
and Mamon & Lokas models on the basis of the radial ve- 
locity dispersion only. Moreover, the differences among these 
models may be partially washed out when considering the 
projected velocity dispersion since this enters also other pa- 
rameters (related to the luminous components) thus leading 
to possible degeneracies that seriously complicate the inter- 
pretation of the data. 

Introducing an anisotropy in the velocity space also 
changes the DF of the model that now becomes a function 
not only of the binding energy £, but also of the modu- 
lus L of the angular momentum vector. The construction 
of general anisotropic DF f(£,L) for spherically symmetric 
potential- density pairs is discussed in detail in Dejonghe 
(1986). However, here, we restrict our attention to the Osip- 
kov - Merrit parameterization in which case the DF depends 
on £ and L only through the parameter : 



Q = £ 



zl OM 



(50) 



It is then possible to compute the DF using the following 
generalization of the Eddington formula (Dejonghe 1986; 
Binney & Tremaine 1987) : 



HQ) 



with 



V8tv 2 Jo 



d 2 pQ d* 



d* 2 y/Q- * 



Pq(t) = 



1 + 



(— y 

\roM J 



p{r) 



(51) 



(52) 



Defining the dimensionless parameter Q = Q/& s and chang- 
ing variable, we get for PoLLS models : 



AO 



P-2 



V8$-27r 2 $-2 

where we have defined : 



X /Oil 



f0M,red 




d pQ dpQ d 
dx 2 dx dx 



(53) 



d* 

dx 



(54) 



with pQ = pq/ps and xq(Q) obtained by solving Eq.(26) 
after replacing £ with Q. Not surprisingly, the reduced DF 
foM,red may not be evaluated analytically, but it is straight- 
forward to compute numerically. The result for the N04 
model is shown in Fig. 13 for three arbitrarily chosen val- 
ues of the anisotropy radius. Comparing with Fig. 6 for the 
isotropic case, we note that the reduced DF is no more ap- 
proximated by a simple power law. The dependence on the 
anisotropy dimensionless radius xom is quite strong for low 
values of Q with the reduced DF being larger for smaller 
xom- On the other hand, there is almost no dependence at 
all on xom for high values of Q since, in this regime, it is 
£ >> L/2r OM so that Q ~ £ and we reduce to the isotropic 
case. Moreover, whatever is the value of vom, the DF is di- 
vergent at the upper end of Q (as in the isotropic case with 
Q replacing £) so that different values of the anisotropy ra- 
dius have a negligible impact. 



7 CONCLUSIONS 

In an attempt to look for underlying similarities among the 
plethora of spherically symmetric galaxy models discussed in 
literature, we have shown that most of the cuspy models may 
be considered as belongings to a general family defined by 
Eq.(2) for the logarithmic slope a(r) of the radial profile. We 
have then restricted our attention to a subclass of models, 
dubbed PoLLS, whose logarithmic slope is a power law in 
r since this is in agreement with the results of the most 
recent and detailed numerical simulations available up to 
now (Navarro et al. 2004). 

PoLLS models are physically meaningful only if the 
slope parameter 7 is positive. The mass density is non sin- 
gular in the centre so that the models are not cuspy, while p 
is an exponentially decreasing function of r so that the total 
mass turns out to be finite. This is an interesting feature 
of the model that makes it different from, e.g., the popular 
NFW model (Navarro et al. 1996, 1997) that has to be trun- 
cated at the virial radius in order to define a finite total mass. 
The mass profile, the rotation curve and the gravitational 
potential of PoLLS models have been evaluated analytically 
and turn out to be expressed in terms of special functions 
that are quite easy to handle (e.g., in numerical codes). In 
particular, while being asymptotically keplerian, the rota- 
tion curve is only slowly declining so that PoLLS models 
may be good candidates for the dark halo of spiral galax- 
ies. Moreover, the position of the peak of the rotation curve 
may be used to constrain the value of the slope parameter 
thus breaking the degeneracy between 7 and the total mass 
that have the same effect on both the rotation curve and the 
velocity dispersion. 

Assuming isotropy in the velocity space, we have nu- 
merically evaluated the velocity dispersion that turns out to 
decline in the outer regions of the system thus offering the 
intriguing possibility that such models may reconcile the 
presence of a dark halo with the apparent lack of dark mat- 
ter suggested by the declining velocity dispersion measured 
in intermediate luminosity elliptical galaxies (Romanovsky 
et al. 2003). Actually, it is not this quantity that has to be 
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compared with the data, but the luminosity weighted line of 
sight velocity dispersion that we have computed assuming 
that the PoLLS model is the only galaxy component. This is, 
of course, a quite unrealistic assumption so that a detailed 
comparison will be performed in a forthcoming paper where 
we will also take care of the surface brightness to constrain 
the modelling of luminous components. 

Being the models spherically symmetric and isotropic, 
the DF only depends on the binding energy £ and may be 
estimated using the Eddington formula. It turns out that 
the DF may be approximated as a single power law over a 
quite large range for £ , but the width of this interval and 
the goodness of the approximation depend critically on the 
values of the slope parameter 7. The density of states has 
also been evaluated for completeness and shows some qual- 
itative features similar to the DF, but with the opposite 
dependence on £ being a decreasing rather than an increas- 
ing function of the binding energy. The product of the DF 
and the density of states gives the differential energy distri- 
bution N that turns out to be a decreasing function with an 
abrupt cutoff at the upper extreme of the range for £ due 
to the finite central density. 

The properties of the PoLLS models have all been inves- 
tigated by assuming both spherical symmetry and isotropy 
in the velocity space. As a first step towards a more re- 
alistic description, we have also considered the effects of 
non vanishing anisotropy parameter /3 CT on the radial ve- 
locity dispersion exploring two different parameterization, 
namely the popular Osipkov- Merrit model (Osipkov 1979; 
Merritt 1985) and the one recently proposed by Mamon & 
Lokas (2004a,b). In both cases, introducing an anisotropy 
in the velocity space leads to an increase of the radial ve- 
locity dispersion with respect to the corresponding isotropic 
models thus originating the well known mass - anisotropy de- 
generacy. Actually, a combined analysis of both the line of 
sight velocity dispersion and the rotation curve could help 
in breaking this degeneracy being v c insensitive to /3 CT . 

We have also estimated the DF for Osipkov - Merrit 
anisotropic PoLLS models finding out that the power law 
approximation for the DF breaks down. As it is expected, 
the DF now depends also on the anisotropy radius tom, but 
this dependence cancels out when considering high values of 
Q since the higher is Q, the more isotropic is the model and 
the nearer is the DF to its divergence at Q ma x- 

PoLLS models are a first attempt towards a different 
approach to galaxy modelling. Rather than starting from an 
expression for the mass density, our initial guess was on the 
logarithmic slope of the radial profile. This quantity gives 
an immediate understanding of the asymptotic behaviours 
of the density profile. Moreover, it seems that a(r) is better 
constrained from numerical simulations. While our atten- 
tion on PoLLS models was mainly motivated by the results 
of Navarro et al. (2004) that refer to dark haloes, modelling 
a galaxy also requires to consider its luminous components. 
These latter may be efficiently described by Zhao models 
(1996, 1997) that represent a generalization of both the pop- 
ular Jaffe (1983) and Hernquist (1990) profiles. Kinematical 
and dynamical data may then be used to constrain the prop- 
erties of the full galaxy modeled using a Zhao model for the 
luminous component and a PoLLS one for the dark halo. 

We would like to conclude with a general comment. As 
the processing speed of modern supercomputers evolves, in- 



troducing new (sophisticated or not) models to fit the re- 
sults of increasingly detailed numerical simulations becomes 
a more and more attractive temptation. Actually, we be- 
lieve that time is come to look for similarities rather than 
discrepancies among the models in order to find out what is 
the right answer to the questions posed by the structure of 
numerically simulated haloes. In our opinion, approaching 
the problem starting from the logarithmic slope may be a 
good road to successfully arrive at the final destination. 
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APPENDIX A: LOGARITHMIC SLOPE FOR 
SOME KNOWN MODELS 

Eq.(2) is a general expression for the logarithmic slope of 
many spherically symmetric models. To give a quantitative 
demonstration of this ansatz, we report in Table 1 the val- 
ues of the parameters (/3, a, b, c, d, 7) for some cuspy models 
discussed in literature. It is worth noting that for all the 
models it is b — d = and c = 1 with only two remark- 
able exceptions. The first is the cuspy halo model (Munoz et 
al. 2001) that is actually used in lensing applications, but 
has not been studied dynamically. The second case is the 
N04 model (Navarro et al. 2004) that belongs to the class of 
PoLLS models we have discussed here. 



APPENDIX B: APPROXIMATED FORMULAE 
FOR THE N04 MODEL 

The basic properties (such as the mass profile, the rotation 
curve and the gravitational potential) of PoLLS models may 
be expressed analytically in terms of special functions, while 
this is not for some other interesting quantities as the ve- 
locity dispersion and the DF. These quantities may be gen- 
erally expressed as the product of a scaling factor that only 
depends on r_2 and M to t and a definite integral that can 
be evaluated numerically for fixed values of the slope pa- 
rameter 7. To save computing time in applications, we give 
here analytical approximated formulae for the N04 model 
for which it is 7 = 0.17. 

Let us first consider the isotropic velocity dispersion 
which may be factorized as 



a 2 (x) 



r-2 ' ' 1 



r(3/ 7 ) 

with la, defined in Eq.(24), to be computed numerically. For 
the N04 model, this latter quantity may be approximated 



log/ CT (a;;7 = 0.17) 



with : 



1 (log a;) 



(Bl) 
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Table 1. Density profile (with x = r/r s and r s a typical scale radius) and values of the parameters (/3, a, b, c, d, 7) for some popular 
spherical galaxy models. In the last column, we also give some useful references. 



Model 


P 





a 


b 


c 


d 


7 


Ref. 


Isothermal 


x- 2 


2 

















Binney & Tremaine 1987 


Jaffc 


x- 2 (l+x)- 2 


2 


2 





1 








Jaffe 1983 


Hcrnquist 


x-^l+x)- 3 


1 


4 





1 








Hernquist 1990 


Denhen 




Id 


4/7D 





1 








Dchncn 1993 


NFW 


x- 1 (l+x)- 2 


1 


3 





1 








Navarro et al. 1996, 1997 


Generalized NFW 


X — @GNFW M + x ^0GN FW— 3 


Pgnfw 


3/PcNFW 





1 








Jiing & Suto 2000 


Cuspy halo 




1CH 





nhcH 





1 





Munoz et al. 2001 


N04 


see Eq.(lO) 


2 














0.17 


Navarro et al. 2004 


RTM 


^"1(1+ ^-3/2 


1 


5/2 





1 








Rasia et al. 2004 



(00,01,02,03) = (7.50404,-2.06006,-0.66146,-0.07771) . 

This approximation works quite well with the error 
AI a /lt rue (being A7 CT = J* r " e - ll u ) which is lower than 
~ 0.5% for < x < 0.1 and even smaller (lower than 
~ 0.08%) in the range (0.1, 10). 

A fundamental quantity to study the dynamical prop- 
erties of a galaxy model is the DF. The reduced DF of the 
isotropic N04 model may be approximated as : 

n— 3 

log f red (£) ~ 1.2052(1 - e)-°' 338 e " e " ( B2 ) 

with e = log£/log£*, = 1.614 and: 

(eo,ei,e 2 ,e 3 ) = (1,0.923,0.0275,0.00476) . 

The error is lower than 4% on / re£ j for —2 < log£ < 0.15 
so that the approximation is indeed quite good. It is worth 
stressing, however, that Eq.(B2) must not be used outside 
the range quoted since leads to serious systematic errors. For 
instance, the approximated reduced DF does not diverge for 
£ — > Emax as it is for the the true reduced DF. 

Finally, let us consider the case of the surface density E 
that we may write as E s x E n with E n = S(y; 7)/<S(l; 7) nu- 
merically evaluated from Eq.(28). It is convenient to use the 
logarithm of E n (that is immediately related to the surface 
brightness) rather E n itself when looking for an approxi- 
mated analytical expression. For the N04 model, we have 
found out that : 

4 

logE n (y; 7 = 0.17)~^ S „(logy) n (B3) 

n=0 

with : 

(so, si, s 2 , s 3 ) = (-0.002, -1.159, -0.369, -0.369, -0.0030) , 

reproduces E n with an error lower than ~ 2.5% over the full 
range y = R/r- 2 G (0,0.1), being AE n /E^ r " e lower than 
~ 0.5% for 0.1 < y < 5. 
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